Hypothalamic transcriptome analysis reveals male-specific differences in molecular pathways related to oxidative phosphorylation between Iberian pig genotypes

The hypothalamus is implicated in controlling feeding and adiposity, besides many other physiological functions, and thus can be of great importance in explaining productive differences between lean and fatty pig breeds. The present study aimed to evaluate the hypothalamic transcriptome of pure Iberian (IBxIB) and Large White x Iberian crossbreds (IBxLW) at 60 days-old, produced in a single maternal environment. Results showed the implication of gender and genotype in the hypothalamic transcriptome, with 51 differentially expressed genes (DEGs) between genotypes and 10 DEGs between genders. Fourteen genotype by sex interactions were found, due to a higher genotype effect on transcriptome found in males. In fact, just 31 DEGs were identified when using only females but 158 using only males. A higher expression of genes related to mitochondrial activity in IBxIB male animals (ND3, ND4, ND5, UQCRC2 and ATP6) was found, which was related to a higher oxidative phosphorylation and greater reactive oxygen species and nitric oxide production. IBxLW male animals showed higher expression of SIRT3 regulator, also related to mitochondrial function. When females were analysed, such differences were not found, since only some differences in genes related to the tricarboxylic acid cycle. Thus, the results indicate a significant effect and interaction of the breed and the sex on the hypothalamic transcriptome at this early age.


Introduction
Currently, intensive pig production relies on lean pig breeds, such as Large White, since these animals have been highly selected for meat production and, therefore, have rapid growth, high

Bioinformatic analysis
FastQC program version 0.11.8 [27] was used to assess the quality of raw sequencing data obtained from the 20 hypothalamic samples. TrimGalore version 0.5.0 [28] was used to qualitative trim data with default settings and to remove the sequencing adaptors and poly A and poly T tails (stringency of 6 bp, -s 6), keeping only paired-end reads where both pairs were longer than 40 bps and that had an optimal Phred Score (Q > 20). Filtered reads were mapped against the pig reference genome Sscrofa11.1 /Ensembl release 94 using HISAT2 version 2.1.0 [29]. Samstools-1. 9 [30] was used to convert the SAM files obtained in the previous step into BAM archives. Read counting and merging was performed with HTSeq-count version 0.11.1 [31].

Differential expression analysis
The files with the counts of reads' number mapped to each gene obtained were analyzed with DESeq2 R package [32] using R 4.0.3. A filtering of the Differential Expressed Genes (DEGs) was obtained under the following two criteria: a Fold Change value (FC) � 1.2 and a Benjami-Holchberg adjusted p value < 0.1. DESeq2 software supports more complex experimental designs in addition to two-groups setups. Using this software, RNA-seq read counts were modelled by generalized linear models, including the genotype and sex effects, the genotype effect within each sex and a full model including genotype, sex and the genotype by sex interaction effects. The comparison of the different DEGs obtained in each analysis was performed using a Venn diagram calculated and draw in http://bioinformatics.psb.ugent.be/webtools/ Venn/. Due to anormal results in one of the individuals (S1 Fig), the DESeq2 analyses were finally performed with the hypothalamus of 19 out of the 20 samples taken (9 from IBxIB piglets; 5 females and 4 males, and 10 from IBxLW piglets, 5 females and 5 males).

Functional interpretation analysis
An enrichment analysis based on the functional annotation of the differentially expressed genes was performed using the Ingenuity Pathway Analysis software (IPA) [33] and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database version 11.0 [34]. In both cases, the different list of DEGs obtained with DESeq2 were uploaded into the software to determine activated pathways, functions, regulator activity or gene correlations [35,36].

Result validation by quantitative PCR (qPCR)
RNA obtained from the 19 animals employed for RNA-seq was used to perform the technical validation of five genes affected by genotype and/or sex (GRPEL2, PDKZ1, SGCA, TBCD and UQCRC2). The genes GAPDH and ACTB were selected for normalization after testing their stability with Genorm software (0.271 < M < 0.363) [37]. Additional information on the selected genes and amplification primers can be found in S2 Table. Primers were designed using Primer Select Software (DNASTAR, Madison, WI, USA) from the available ENSEMBL sequences. Quantitave PCR (qPCR) was performed as previously described [37].
The technical validation was performed by studying the Pearson correlation between the expression values obtained from RNA-Seq data and the normalized gene expression data obtained by RT-qPCR. To validate the global RNA-Seq results, the concordance correlation coefficient (CCC) [37] was calculated between the FC values estimated from RNA-Seq and qPCR expression measures for the five genes. Statistical analysis was performed with a linear model fitting genotype, sex and the genotype by sex interaction as fixed effects, and litter as random effect. All the analyses were performed using MIXED procedure of SAS 9.4 (SAS Inst. Inc., Cary, NC, USA). The model used was: Being y the qPCR expression result for each gene; i = 2 levels (pure Iberian pig and Iberian x Large White crossbred), j = 2 levels (female and male); k = 10 levels (number of litters).

Results and discussion
In a previous study, we assessed the phenotypic data of the animals used in the present research [24] corresponding to the whole generated population. Firstly, from birth to weaning, body measures and weight were taken to ascertain development during lactation. After selecting 67 piglets, plasma metabolic and pro/antioxidant indexes and fatty acid composition of various tissues, as well as postnatal average daily weight gain and growth was assessed at different time-points.
In brief, IBxLW piglets were significantly heavier and larger than their IBxIB counterparts, with differences in body-weight and body-length being more prominent between males than between females. IBxLW showed higher meat content, whereas IBxIB animals had higher adipose tissue content in terms of both, subcutaneous and intramuscular fat, with higher MUFA and SFA content and less PUFA. Organs were also heavier in IBxLW piglets, but some ratios, such as brain/head ratio, were higher in IBxIB animals. At 60 days old, antioxidant capacity and lipid peroxidation levels were also different between genotypes, with the first being affected by a genotype by sex interaction (both IBxIB females and IBxLW males had greater values than their same-genotype counterparts) and the latter being higher in the IBxLW group. Fructosamine was also affected by a genotype by sex interaction, and lactate levels were significantly higher in IBxLW piglets. Those results were also found in the 19 piglets selected in the present study (Table 1 and S3 Table), demonstrating the representative nature of the sample.

Mapping and annotation
From the hypothalamic transcriptome of 20 young piglets, an average of 47.67 million reads were obtained per sample, ranging from 34.95 to 60.50. Read length was 76 bases, and quality (Phred) Score was approximately 40. GC content ranged from 49 to 56%. HISAT2 mapping resulted in a mean alignment rate of 97%, which is in accordance with previous studies by our group in muscle [38] and in adipose tissue [39].

Genotype effect: Differential expression and functional analysis
Firstly, the differences in the gene expression between the hypothalamic samples of IBxIB (excepting the one excluded for being an outlier; 9 piglets) and the total number of samples from IBxLW piglets (10 animals) were studied (both genders together in each genotype, so these results are independent of sex). A total of 51 differentially expressed genes (DEGs) were found. Eighteen were overexpressed in IBxIB animals (log2FC < -0.26) and 33 were overexpressed in IBxLW (log2FC > 0.26). Table 2 contains the top ten differentially expressed genes (5 overexpressed in IBxIB and 5 in IBxLW piglets), whereas the full detailed list of DEGs between both groups can be found in S4 Table. The gene showing the highest overexpression in the IBxIB group was WFIKKN1 (log2FC = -1.88; B-H adjusted p value = 0.09). WFIKKN1 encodes a secreted multidomain protein with implications in inhibition of proteases [40] acting as a regulator and as a potent antagonist for different growth and differentiation [41][42][43]. Therefore, its upregulation in IBxIB animals could be related to the lower development found in these piglets when compared to their IBxLW littermates.
Regarding IBxLW animals, KY was the most over-expressed gene (log2FC = 2.53, B-H adjusted p value = 0.099). It is implicated in the function, maturation and stabilization of the neuromuscular junction and, possibly, in the normal muscle growth [44]. Therefore, the finding of this gene being upregulated in IBxLW when compared to IBxIB pigs may be related to the lower growth potential of the Iberian breed [4], although further research would be needed to clarify the exact function of this gene at the hypothalamus level. Interestingly, leptin receptor (LEPR) gene was not differentially expressed between genotypes. Different expression levels between lean and Iberian pigs, with the latter having significantly lower expression levels, is a usual finding in studies comparing pure Iberian and crossbreds in tissues such as muscle [23] and liver [45] at fetal stages, whereas higher expression levels of LEP can be found in Iberian pigs compared to other breeds [46]. Ovilo et al [17] also found differences in the expression of this gene at the hypothalamic level which were associated with a missense polymorphism in the LEPR gene (LEPRc.1987T/C). Allele LEPRc. 1987T, fixed in the Iberian breed is associated with a lower hypothalamic expression of the leptin receptor. This polymorphism was used for genotype confirmation in the animals studied in the present work, with Iberians being TT and crossbreds being CT, thus differential expression of the gene was expected. A possible explanation of the result obtained would be that the age could be an important factor influencing the results. However, further research in this regard should be performed to fully understand this outcome.
When the DEGs obtained in the genotype comparison employing all piglets (S4 Table) was used for functional interpretation, results with STRING software were only related to Ribonuclease CAF1 and to the CCR4-NOT transcription complex subunit 7/8/Pop2 implicated in the 3' to 5' mRNA deadenylation [47, 48], a process related to synaptic changes and plasticity and inflammatory processes [49], which is consistent with our results. Using IPA, limited results were obtained, with no differences in canonical pathways, the upstream analysis nor regulators. There were differences only in specific functions, such as Nervous system development and function, in which 5 molecules were implicated (SLC6A9, ADCYAPI, ITPRI, MAPK7 and LSS), with p-values ranging from 4.36x10 -2 to 1.86x10 -3 .

Sex effect: Differential expression and functional analysis
When differences between sexes (by the comparison of IBxIB and IBxLW females vs. IBxIB and IBxLW males) were analysed, 10 genes were differentially expressed (Table 3). Only RBBP7 gene was annotated, the remaining being novel genes. Two genes were overexpressed in females (including RBBP7), the rest of them being overexpressed in males. No results were obtained in the functional analysis.

PLOS ONE
Differences in the hypothalamic transcriptome between Iberian pig genotypes RBBP7 is a nuclear protein that belongs to a highly conserved subfamily of WD-repeat proteins that can be found in many histone acetyltransferase complexes [50]. Its principal function is to regulate cell differentiation and proliferation by binding to retinoblastoma protein [51]. At the central nervous system level, low expression levels of this gene have been related to alterations in brain including Alzheimer's disease, neuritic plaque density and Braak Staging [52,53]. In the present study, the higher expression found in females compared to males could be related to the general higher survival rate found in female piglet when compared to male counterparts [54-56] that could be explained by a higher or earlier neural development in this sex [57].

Genotype by sex interaction effect: Differential expression and functional analysis
Regarding the genotype by sex interaction, there were 14 qualitative significant interactions ( Table 4), indicating that, even if, as stated before, sex alone had limited effect on our results, it was important when considered with the genotype effect.
The highest fold change was observed for GRPEL2. This molecule is highly influenced by oxidative stress, being a "sensor", that is capable of augment the activity of mtHsp70 in cases of high oxidative stress levels to prevent the misfolding of mitochondrial imported proteins [58, 59]. The higher expression observed in IBxIB males and IBxLW females is in accordance with the results found in the malondialdehyde (MDA) levels when analysing the phenotype of these animals [24], since IBxIB males and IBxLW females had higher MDA levels than their same

Genotype effect within sex: Differential expression and functional analysis
Phenotype analysis of the animals showed a strong sex effect, with differences between IBxIB and IBxLW males being higher than when comparing females [24]. Furthermore, a great number of genotype by sex interactions were found when performing the differential expression analysis with DESeq2. Thus, to fully understand the processes occurring in males and in females in the hypothalamus that could explain the phenotypic differences, the 19 piglets were divided according to their sex, and the genotype effect was studied in each sex separately.
When only females were studied, 31 DEGs were found according to genotype. Thirteen DEGs were overexpressed in IBxIB, and 18 in IBxLW female piglets (Table 5 and S5 Table).
FAM160B2 (log2FC = -1.12; B-H adjusted p value = 0.082) and USP16 (log2FC = 1.82; B-H adjusted p value = 0.025) were the genes with the highest Fold Change in IBxIB and IBxLW females, respectively. FAM160B2, also known as Retinoic acid induced 16 (RAI16), is a gene possibly engaged in the maintenance of the intestinal barrier and in the immunoprotective inflammation due to its ability to activate the MAPS/ERK and TGFB signalling pathways [61]. It has been studied at the hepatic tissue [62] and in prostate cells [63]. However, at the central nervous system (CNS) only its analogous FAM160B1 has been studied, with results indicating that deficits of this gene are related to severe microencephaly and intellectual disability [64]. Thus, further studies for FAM160B2 at CNS level are necessary to ascertain possible implications in its development and, therefore, survival of the animal.
On the other hand, USP16 is implicated in T-cell mediated immunity [65], and also in the LDL receptor stability and activity [66]. The latter function could have relevant implications in the piglets of the present study. The Iberian pig, with its natural leptin resistance, shows similar characteristics to those of obese and diabetic humans, including high levels of LDL-c in blood in the later age [8,9]. Thus, the higher expression of this gene in IBxLW females could be related to the overall lower levels of LDL-c in this group, even if at 60 days-old no significant differences were found [24]. When functional analysis was performed with the set of DEGs detected in only females, results in STRING showed a reactome pathway related to the tricarboxilic acid (TCA) cycle and respiratory electron transport, with two genes implicated (CHD8 and SLA-DQB1; overexpressed in IBxIB females). Since both genes implicated in this pathway were overexpressed in IBxIB females, this result could mean higher oxidation rates in pure Iberian animals than crossbreds and, therefore, an impairment of the modulation of energy balance at the hypothalamic level [67], with possibly higher production of reactive oxygen species (ROS) during the oxidative phosphorylation reactions, since both processes are linked [68].
The genotype effect was also evaluated within males. Accordingly to their higher differences in phenotype characteristics [24], IBxIB and IBxLW males showed a greater number of DEGs than females (158 DEGs, 120 overexpressed in IBxIB and 38 in IBxLW males). Excluding a novel gene, CDADC1 and TBCD were the two genes with the highest log2FC (-1.04 and 1.32, respectively; Table 6 and S6 Table).
CDADC1 (or NYD-SP15) is a protein coding gene, highly overexpressed in IBxIB males. It has a potential inhibitory effect on cell growth [69], and also has important implications in spermatogenesis and testicular development [70]. This result is in accordance to CRISP1 (also related to male fertility, [71]) overexpression observed in pure animals when the total number of animals (IBxIB vs. IBxLW) were analysed. Thus, both findings could be related to an earlier onset of puberty in IBxIB animals. However, it should be highlighted that these piglets were castrated during the first week of life, so further research is needed to confirm this hypothesis.
TBCD, overexpressed in IBxLW pigs, is pivotal in the correct folding of microtubulins [72]. Therefore, it plays a major role in the adequate function and development of the body, especially the nervous system, since mutations in this gene have been related to different neuropathologies [73,74]. This outcome agrees with results obtained at the phenotype level. Although IBxLW females also overexpressed this gene, they did it in a lesser extent than males. As explained before, differences between males in body-weight and size were greater than between females [24], an outcome in which TBCD gene could be implicated. This would be in accordance with previous results in Iberian pigs in which a higher growth rate was seen in early development in females when compared to males in IUGR conditions [56]. The higher growth of females in the early development compared to males has also been described in humans [75] and rats [76,77], where a different feeding behaviour between females and males

PLOS ONE
Differences in the hypothalamic transcriptome between Iberian pig genotypes controlled by the hypothalamus was hypothesized. Even though the animals selected in the present study were not IUGRs animals, it is plausible that, being a rustic breed, Iberian females may be better prepared for survival, given their importance in the species continuity. This would explain the limited differences observed in females at the transcriptome level in comparison to males, with a lower number of DEGs or, in the case of some essential genes for development as TBCD, with the difference in expression being less significant. The list of male-specific DEGs was analysed to study its functional implications. In this case, processes related to mitochondrial respiration, electron transport and oxidative phosphorylation showed a high relevance, with a great number of DEGs implicated in such biological functions. STRING software allowed the clustering of functionally related DEGs involved in these processes (Fig 1).
Results obtained with the software IPA also showed significant differences in functions such as Cellular Assembly and Organization, Cell-To-Cell Signalling and Interactions, Lipid Metabolism or Organismal Development (number of molecules ranging from 16 to 40, p value ranging from 6.87x10 -3 to 3.82x10 -5 ). Ten regulators showed significant z-scores ( Table 7). Five of them were activated in the IBxIB genotype (negative z-score), whereas the other five were activated in IBxLW animals (positive z-score).
When analysing IBxIB vs. IBxLW males, mitochondrial gene regulators such as LONP1, ALKBH1, NSUN3 or DAP3 showed negative significant z-scores (activated in IBxIB males), and the canonical pathway Oxidative Phosphorylation showed a significant z-score (z-score = -2.121; overrepresented in IBxIB animals), implying a higher activation of this pathway in pure Iberian males. In fact, the genes implicated in the canonical pathway were the same highlighted by the STRING software. All mitochondrial genes except NDUFS2 were overexpressed in IBxIB animals (Table 8). In fact, the hypothalamic mitochondrial gene regulation found in this research was similar to the one found in the muscle in previous studies, since these results are in accordance with recent work reporting gene expression differences between lean and fatty Serbian pig breeds analysing the muscle transcriptome of adult animals [37], in which oxidative phosphorylation was the main process affected by genotype, being activated in the muscle of the fat Mangalitsa breed.
Oxidative phosphorylation (OXPHOS) is the process in which ATP is formed because of the transfer of electrons from NADH or FADH2 to O2 by a series of electron carriers (complexes). It constitutes the major source of ATP in aerobic organisms, being performed in the inner layer of the mitochondria [78]. Fig 2 illustrates this process, with genes overexpressed in IBxIB or IBxLW males coloured in green or red, respectively. Concretely, ND3, ND4 and ND5 (overexpressed in IBxIB pigs) and NDUFS2 (overexpressed in IBxLW pigs) are implicated in the correct function of the complex I [79], which is the major entry point for electrons. Thus, Table 7. Regulators predicted by IPA to be activated in 60 days-old pure Iberians (IBxIB; negative z-score) and in Iberian x Large White crossbred pigs (IBxLW; positive z-score). These regulators potentially control the differences in the hypothalamic transcriptome observed in males.

Molecule Type Predicted Activation State
Activation zscore p-value of overlap

PLOS ONE
Differences in the hypothalamic transcriptome between Iberian pig genotypes it is considered the rate-limiting step in mitochondrial respiration, playing an essential role in energy metabolism [80]. UQCRC2 is a core protein of the complex III [81], which is overexpressed in IBxIB males. In this complex, electrons, after passing through the coenzyme Q, are transferred from one subunit to other to produce ATP [82]. Afterwards, electrons are transferred to complex IV, being the DEGs COX2 and COX6 (more expressed in IBxIB males) part of this complex. The electrons are transferred to oxygen and hydrogen (with high levels of energy yielded) [83], whereas protons are pumped across the membrane [84]. ATP6 is a DEG also with higher expression levels in pure males than crossbreds involved in the complex V, the last step of the respiratory chain, which transforms ADP to ATP in the mitochondrial matrix using the energy provided in the previous steps of the process [85]. OXPHOS has been related to obesity and type 2 diabetes, being considered a key factor in the development of these disorders [86][87][88]. The excessive intake of nutrients (or leptin resistance states) makes the mitochondria overloaded with nutrients such as glucose and fatty acids, which increases the production of Acetil-CoA, NADH and, therefore, the proportions of electrons that enters the mitochoncrial intermembrane space [89]. This process is a major productor of reactive oxygen species (ROS; reviewed in [90]). In fact, another canonical pathway affected by genotype in males was Production of Nitric Oxide and Reactive Oxygen Species in Macrophages (enriched in IBxIB males when compared with same-sex piglets from the IBxLW group; z-core = -1.000). Both, nitric oxyde (NO) and ROS lead to oxydative stress [91], which is related to inflammation, adiposity, metabolic syndrome and insulin resistance [92,93]. Specifically, the production of NO is closely associated with the hypothalamic inflammation in obese states, being mediated by macrophages [94]. These inmune cells play a major role in the diabetes progression, since they are also major producers of ROS [95]. ROS also mediates the higher differentiation of macrophages to M1 type, which are proinflammatory, and implicated in the chronic inflammation found in obese and diabetic individuals [96]. These results obtained at the hypothalamic level in the present study are in agreement with the outcomes found when comparing the muscle transcriptome of Mangalitsa and Moravka Pigs [37], in which genes related to OXPHOS and mytochondrial dysfunction were also upregulated in the fatty breed (Mangalitsa) when compared to lean ones (Moravka), even if the age, breeds and tissues are different between both studies. Thus, results seemed to be more probably related to the obesity and lean-type of each breed than the genotype itself. This fact could implicate a relation of the increased phosphorylation with obesity and diabetoid states, as seen in humans [97,98].
Results from IPA also showed potential differences in important metabolism regulators such as STAT3 (z-score = -1.000; its causal network is represented in Fig 3). The STAT is a family of transcription factors with important functions in signal regulation [99][100][101]. Concretely, STAT3 has multiple isoforms, and regulatory effects [102]. This regulator is being increasingly studied in the hypothalamus because of its important relation with leptin levels and, therefore, leptin resistance (as reviewed by Liu et al [103]). Thus, in response to leptin, a transcriptional cascade is activated, including STAT3, which has important effects on energy balance [104]. At the hypothalamus level, the stimulation by leptin makes STAT3 bind proopiomelanocortin (POMC) [105]. POMC is a propeptide expressed in the hypothalamus which produces in this location the rising of different substance, such as α and β-melanocyte stimulating hormones. Those hormones have anorexigenic functions, also enhancing insuling sentitivity and stimulating energy expenditure [106]. Thus, a lower POMC expression leads to severe hyperphagia and obesity in humans, mice and pig [107][108][109]. Obesity produces hypothalamic inflammation and higher stress, impairing, in turn, POMC function [110][111][112]. Accordingly to the greater ingestion capacity of the Iberian pig compared to commercial breeds [113], the POMC gene was overexpressed in IBxLW in comparison to IBxIB males ( Table 6, red node in Fig 3). This result is in agreement with the leptin resistance syndrome, characteristic of Iberian pigs, because a reduced leptin signalling leads to a lower expression of downstream anorexigenic genes such as POMC, although, as explained before, no expression difference was observed at the LEPR gene level. Also, as represented in the figure, STAT3 has important implications in the mitochondrial gene expression and oxidative metabolism [114], which agrees with results regarding differences in oxidative stress between IBxIB and IBxLW stated before. All these outcomes are related to obesity and inflammatory response, being activated in IBxIB animals, as previously demonstrated [8,46].
The canonical pathway Sirtuin Signaling Pathway was potentially enriched in IBxLW male piglets when compared with IBxIB males (z-score = 0.447). Sirtuins are a protein family highly conserved from bacteria to mammals and implicated in the regulation of many biological processes, such as stress, metabolism, development or longevity [115]. Concretely, SIRT3, a mitochondrial sirtuin [116], was detected by IPA as an upstream regulator activated in IBxLW animals (z-score = 2.236). It has been suggested that SIRT3 has important implications in mitochondrial activity regulation, especially by deacetylation of different molecules [117]. This enzyme rises during fasting or caloric restriction [118][119][120][121], and declines as response to highfat diets or in insulin-resistant states [119,122,123]. Since diet in both groups was the same throughout the experiment, it is possible that the difference in SIRT3 activation is associated to the natural insulin-resistant state of the Iberian pig [8], which would explain why it is more activated in crossbreds than in pure males.
All in all, results found when only comparing males were mostly correlated with mitochondrial dysfunction (leading to higher OXPHOS and, therefore, with oxidative stress and inflammation mediated by proinflammatory macrophages), and with differences in metabolism regulated by STAT3 and sirtuins. However, further research is required to fully clarify the implications of these findings in the differences between lean and obese individuals.

Comparative analysis of results among datasets
A comparative analysis was performed between results obtained in the studies of the genotype effect, the genotype effect within each sex and the interaction effect to fully understand them, which is summarized in Fig 4. Most genes were unique in each dataset, which further states the idea of the relevant differences found in each sex in the present study. However, the comparison of all IBxIB vs all IBxLW pigs shared 8 genes with both the analyses using only females and only males (TBCD, TMEM138, CYR61, SLC6A9, DUT, TMEM9B, TACC3 and SGCA). These 8 genes are those for which the genotype effect is clearest and independent of sex. The overall genotype analysis with all animals also shared 9 genes with the comparison IBxIB vs. IBxLW females (PRR19, PGGHG, CHD8, USP16, GYPCÂ, SEPTIN1, TARBP1, CNOT7 and FAM160B2) and other 9 with the analysis comparing only males (WDR45B, TMEM120B,  LRRC45, CCDC84, CPSF7, PDZK1, ADCYAP1, KRR1 and SMAP1). Moreover, the datasets from the comparison of IBxIB vs. IBxLW males and the interaction also had 8 genes in common (LIN37, RBM12B, PLPP1, CELSR1, PSAP, DYNC1H1, GRPEL2 and TMEM165), for which the evidence of a male-specific effect is stronger. Some of them had, as stated before, important implications in explaining phenotypic differences between IBxIB and IBxLW animals, like SGCA, GRPEL2 or TBCD, which may explain the coincidences among analyses.

Result validation with qPCR
To validate the RNA-seq results, the relative expression of a selected group of DEGs was assessed (S7 Table). Due to the particularities of the results obtained, with substantial effect of sex in the differences of the transcriptome between genotypes, the genes selected for validation were differentially expressed in one of the following analyses: all animals, females, or males. Therefore, the RNA-seq and real time qPCR correlation was calculated separately for each subsampling group. As expected, in the group where the gene showed the higher difference of expression, the correlation showed high values, whereas poorer results were obtained when the difference of expression in the RNA-seq analyses was lower or non-significant.
As found in other studies [38,124], differences tend to be more notable with the RNA-seq approach than with the Real-Time qPCR. However, the concordance correlation coefficient (CCC) showed an acceptable correlation between techniques (CCC in all samples = 0.73; CCC in females = 0.66; CCC in males = 0.60).

Conclusions
Differences between 60 days-old pure Iberian and Iberian x Large White crossbreds were found at the hypothalamic transcriptome level, even when bred under the same prenatal and early postnatal environments. In accordance with the phenotype data, genotype effect on transcriptome was greater in males than in females. Given the importance of the hypothalamus in growth patterns and adiposity, these findings help to better understand the metabolic and physiological basis of phenotypic differences between lean and fatty pigs, which are important from a productive and biomedical point of view.